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CN . Abstract 

We present the convergence rates and the explicit error bounds of Hill's method, which is a numerical 
method for computing the spectra of ordinary differential operators with periodic coefficients. This method 
approximates the operator by a finite dimensional matrix. On the assumption that the operator is self- 
adjoint, it is shown that, under some conditions, we can obtain the convergence rates of eigenvalues with 
respect to the dimension and the explicit error bounds. Numerical examples demonstrate that we can verify 
these conditions using Gershgorin's theorem for some real problems. Main theorems are proved using the 
Dunford integrals which project an eigenvector to the corresponding eigenspace. 
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1 Introduction 



This paper considers Hill's method [5], which is a numerical method for computing spectra of ordinary differential 
J> , operators S p with periodic coefficients: 
OO 



, where x € R, and the functions fj (j = 0, 1, . . . ,p — 1) are C°° and satisfy fj(x + L) = fj(x) for some real 
' positive constant L. The eigenvalue problem for S p is described as S p (f> = \<f), where A is an eigenvalue and </> 
is an eigenvector. This eigenvalue problem often appears in physical problems such as linear stability analysis 
of periodic solutions of nonlinear wave equations [7J. The set of eigenvalues A is included in the spectrum 
a(S p ). Accordingly, computation of &(S P ) is important in both theoretical and practical points of view. It was 
reported in [7] that Hill's method can produce very good computed results of cr(S p ) for some problems. Note 
, that implementation of Hill's method is straightforward. 

The ideas of Hill's method are to employ the Floquet-Bloch decomposition of the spectrum cr(S p ) and 
to approximate the eigenvector <f> by a finite Fourier series. This Fourier series approximation generates an 
eigenvalue problem of finite dimension D corresponding to the operator S p . The aims of this paper are to show 
the convergence rates of the approximate eigenvalues with respect to the dimension D, and to explicitly obtain 
the error bounds of them. 

On the convergence property of Hill's method, Curtis and Deconinck [J proved that the exact eigenvalues 
exist near the computed ones for general cases of S p , and showed that all exact eigenvalues can be approximated 
by the computed ones for the case of self-adjoint S p . Also, they obtained the convergence rate for the case of 
self-adjoint S p and the constant coefficients fx, ... , f p -i- Johnson and Zumbrun |10j investigated Hill's method 
using the Evans function, of which the roots correspond to eigenvalues, for general cases of S p . They showed that 
the approximate eigenvalues converge to the exact ones, but did not get the convergence rate of them. Vainikko 
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jll)|14j examined an approximation method for the eigenvalue problem S p <j> = \<j) in an abstract framework, 
and obtained the convergence rate, which is based on the resolvent norm convergence of approximate operators, 
for general cases of S p . We can apply Vainikko's results to Hill's method. But it is difficult to directly evaluate 
the value of the convergence rate. It should be noted that, although all of these convergence rates yield some 
error bounds of approximate eigenvalues with unknown coefficients, any explicit error bounds have not been 
shown. 

In this paper, we give a priori estimates of the convergence rate and a posteriori explicit error bounds of 
Hill's method for self-adjoint S p with the two cases of coefficient functions fj in namely fj G C°° and fj 

being analytic on some strip region containing the real line. These classes of fj are more general than [3]. The 
key ideas of these estimations are to project eigenvectors using the Dunford integrals (|5.1[) and (|5.20[) . and to 
specify the disks around the exact eigenvalues using Gershgorin's theorem (Theorem 13. 3p . 

This paper is organized as follows. Section [2] describes Hill's method. In Section [3j we give the convergence 
rates and the explicit error bounds of Hill's method for self-adjoint differential operators. Section [4] presents 
numerical examples which supports our results. Section [5] summarizes the proofs of theorems in Section [3] 
Section |5] concludes this paper. 



2 Hill's method 

This section describes Hill's method [5] which is a numerical method for computing the spectum (j(S p ) of the 
ordinary differential operator S p with periodic coefficients fj defined by (jl.ll) . This operator S p can be regarded 
as an operator S p : _ff p (R) — > L2(R), where L2CR) is the Lebesgue space of square integrablc functions on R 
and H P (IV) C L2(R) is the Sobolcv space of functions whose derivatives up to p-th order arc square integrablc. 
Hill's method approximates the elements of <j(S p ) by the following two steps. 



Step 1: Floquet-Bloch decomposition. In order to apply the Floquet theory, we introduce a new operator 
Sg : HP([0,L]) per -> L 2 {[0,L]) peI defined by 

S£ = e-^ x S p e^ x , (2.1) 

where p G [0, 2n/L), and L2QO, L]) pcl and H p ([0, L]) pcl - are the Lebesgue space and the Sobolev space of periodic 
functions on [0,L], respectively. Note that S p is defined on R, whereas is on [0,L]. If S p is self-adjoint, then 
Sp 1 is also self-adjoint. Moreover, Sp 1 is explicitly written as 

for some periodic functions /o, . . . , f p -\ with period L. 

It is known that the spectrum cr(Sg) consists of only the eigenvalues of Sp 1 , and the Floquet theory yields 

v(S p )= |J a(Sp. (2.3) 

(lE[0, 2n/L) 

This decomposition (|2.3[) is called the Floquet-Bloch decomposition. Accordingly, it suffices to consider the 
eigenvalue problem S^cf) = \<fi for fj, g [0, 2%/L), where A G C is an eigenvalue and <fi G i2([0, L]) pei is an 
eigenvector. In addition, it should be noted that <j(S^) is a discrete set of the eigenvalues of S£ without ac- 
cumulation points. This discreteness follows from the compactness of the resolvent of Sp 1 provided p(S^) ^ 0, 
where p(S£) = C \ er(S^). Sec e.g. Lemma 2 in [6l Chapter XIX Section 2], Lemma 3 in [6l Chapter XIX 
Section 3] and Lemma 16 in [5j Chapter XIII Section 2]. 



Step 2: Fourier series approximation. Since <j> G L2([0, L]) per has the Fourier series expansion: 



1 -r-^ - ( .2imx\ 



(2.4) 



2 



can be approximated by the truncation of this series: 

(P N <f>)(x) = -j= 0ncxp^-i-^j. (2.5) 



= -N 

This truncation reduces the eigenvalue problem S£4> — \<j) to a finite dimensional problem. More precisely, the 
problem 

S^ N <t>N = \ N <t>N with S£ N = P N S£P N (2.6) 

gives approximate eigenvalues Ajv's for the original problem S£4> — X(f>. 

Since the problem (|2.6[) is equivalent to a matrix eigenvalue problem, we can obtain the approximate values 
of the eigenvalues of using some standard numerical method. In the following sections, let cr(S£ N ) denote 
the set of eigenvalues of S£ N . 

3 Convergence Rates and Error Bounds of Hill's Method 

In this section, we present theorems about a priori estimates of the convergence rate and a posteriori explicit 
error bounds of Hill's method for the eigenvalue problem S£<f> = A<fi with self-adjoint S£ on some assumptions. 
Their proofs are given in Section [5J In what follows, we use the notations defined in Section [5J 

3.1 Assumptions 

First, we assume the self-adjointness of the operator S£. 

Assumption 1. The operator S£ is self-adjoint, i.e. (S£)* = Sg, where * represents the adjoint of operators. 

On Assumption [2 N in (|2.6[) is also self-adjoint. 

Next, we specify the smoothness of the coefficient functions fj's of in (|2.2[) . 

Assumption 2a. For the operator S£, the coefficients /o,/i, • • • , f P -i are C°° . 

In addition, as a special case of Assumption I2a[ the following assumption is prepared. 

Assumption 2b. For the operator S£ , the coefficients /o, /i, . . . , / P -i are analytic on a complex domain 

V d = {z<EC\ |Imz| < d} (3.1) 

for some d > 0. 

Assumption I2bl is often satisfied in real problems. 

If S£ N approximates Sg very well, we may expect that approximate eigenvalues Aat's for large enough N 
should be included in the neighborhood B\(r\) of the corresponding exact eigenvalue A, where 

B ( (r) = {z€C\\z-C\<r}. (3.2) 

Then we prepare the following two assumptions. 

Assumption 3a. For A £ o~(S£), there exists a positive real number r\ such that the following holds true: for 
some sequence {An | An € &(Sp N )} and some positive integer No, if N > No then 

|A-A N |<^ (3.3) 

and 

2r x < min {dist(A, <r(Sg) \ {A}), dist(A, <t(S» n ) \ {\ N })} . (3.4) 
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Figure 1: Illustration of the conditions (|3.3[) and p.4[) in Assumption I3al This figure depicts locations of an 
exact eigenvalue A £ o~(Sg), the corresponding approximate eigenvalue Xn € f(S'p w ) and some other eigenvalues 
on the complex plane C. The eigenvalue Xn is located inside the circle with center A and radius r\/2 (shown 
by the smaller dotted circle), and the other eigenvalues than A and Xn are outside the circle with center A and 
radius 2r\ (shown by the larger dotted circle). The directed contour r\(Sg) is used to prove the theorems in 
Sections 13.21 and 13.31 



Figure Q] illustrates these conditions (|3.3[) and (|3.4[) . When a single sequence {Xn} for some A £ o-(Sg) 
is detected, Assumption l3al is effective for estimation of the convergence rate and the error bound of {Aat}. 
In general, there exist more than one sequence of approximate eigenvalues for some A £ o~(S£). Then we set 
another assumption as follows: 

Assumption 3b. Let k be an integer with k > 2. For X £ a(S^), there exists a positive real number r\ such 
that the following holds true: for some sequences {Xn,% | Xn,% £ &(Sp N )} {i = 1,2, ...fc) and some positive 
integer No, if N > Nq then 

|A-Ajv,i|<Y (i = M,...,fc) (3.5) 

and 

2r x < min |dist(A, <r(S£) \ {A}), dist ^A, a(S^ N ) \ ( LH J } ' ( 3 ' 6 ) 

Remark 3.1. Assumptions [3a\ and\3b\ mean that, for sufficiently large N, the number of eigenvalues of S p l N 
approximating X £ o~(S£) is at most finite, and the other elements of o~(S£) and o~(Sp N ) are relatively far from 
X. The finiteness of the approximate eigenvalues is based on the fact that the eigenspace of for any X is finite 
dimensional. 

When Assumption I3al or I3bl arc satisfied, there exist some convergent subsequences in [J N o~{S^ N ). Thus 
these assumptions are closely related to the convergence property of the eigenvalues of Sr n - On this conver- 
gence property, the following two theorems are known [J] [T^] . Theorem 13.11 gives so-called "no-spurious- mode 
condition" . 

Theorem 3.1 ([5J Theorem 9]). Let D C C be a compact set and {Aat} be a sequence in 2) with Xn £ o~(Sp N ). 
Then for any e > there exists a positive integer M such that the following holds true: for any N with N > M , 
Xn is contained in the e -neighborhood of some X £ 23 f~l o~(S£). 

Theorem 3.2 ([H Corollary of Lemma 10 and Theorem 16]). Assume that S£ is self-adjoint. Then for any 
X £ o-(S^) there exists some subsequence {Ajv^} of {Xn} with Xn £ °~{Sp N ) such that Xnj —> X. 

In Section 13.31 a sufficient condition for (|3.3p and (|3.4p in Assumption [3a| is given by Theorem 13.61 We can 
verify this sufficient condition using Theorems 13. 1 1 13.21 and the following theorem for a concrete problem as 
shown in Section [U 
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Theorem 3.3 (Gershgorin's Theorem [U]). For any n x n matrix V = (wy) with Vij G C, aZZ eigenvalues of 
V are contained in 1J" =1 Ci, where C\, . . . , C n are tfofcs defined as 



d = < jz G C 



Moreover, each connected component o/lJ™ =1 Ci contains as many eigenvalues ofV as the disks composing it. 

This theorem is used to estimate the radius r \ in the conditions (|3.3[) and (|3.4p , whereas Theorems 13.11 and 13.21 
guarantee that the sequence of the approximate eigenvalues converges to an exact eigenvalue. 

3.2 A Priori Estimates of Convergence Rates 

When Assumptions [U 2 and 3 are satisfied, we can obtain a priori estimates of the convergence rates as follows. 

Theorem 3.4. Let Assumptions^ and [3a\ be satisfied. Then, when Assumvtion \2a\ is satisfied, namely fj being 
C°° , for a sufficiently large N and any positive integer q, there exists a positive constant C\ such that 

|A-Ajv| < CiiV"*, (3.8) 

where C\ depends only on p, /o, . . . , f P -i, A, r\, q and an eigenvector <fi of S£ corresponding to A. In addition, 
when Assumption\2$ is satisfied, namely fj being analytic, for a sufficiently large N and any e with < e < d, 
there exists a positive constant C 2 such that 



I A - A* I < C 2 exp (-^pl N ^j 



(3.9) 



where C 2 depends only on p, fo, . . . , f v -\, A, T\, d, e and <f>. 

Theorem 3.5. Let Assumvtions\l\ andUFb] be satisfied. Then, when Assumvtion \2a\ is satisfied, namely fj being 
C°° , for a sufficiently large N , some choice of an index iff G {1, . . . , k} and any positive integer q, there exists 
a positive constant C\ such that 

|A-A W , JJV | <dN-", (3.10) 

where C\ depends only on p, /o, . . . , f p —i, A, r\, q and an eigenvector <f> of S£ corresponding to A. In addition, 
when Assumption\2^ is satisfied, namely fj being analytic, for a sufficiently large N, some choice of an index 
in G {1, . . . , k) and any e with < e < d, there exists a positive constant C2 such that 

|A - X N , N I < C 2 A^ +1 exp (- Md k ~ £) N) , (3.11) 

where C 2 depends only on p, fo, . . . , f p -i, X, r\, d, e and <j>. 

3.3 A Posteriori Error Bounds 

When Assumptions [1] and [2a| arc satisfied, we can explicitly obtain a posteriori estimates of the error bounds 
as follows. Here, for simplicity, only the case of Assumption [3a| is considered. 

Theorem 3.6. Let Assumptions^ and\2a[ be satisfied. Assume that there exists Q G C and r > such that, for 
any N, 

X N G B ( (r) and B c (9r) n (*(S» N ) \ {Ajy}) = 0, (3.12) 
where Xn G cr(Sp N ) and B^(r) is defined by (|3.2[) . Then there uniquely exists X G cr(S^) (1 B^(r) satisfying 

P I l+N \ 

'E E E \(fj)U4>N)l-m\+(2N + l) \(fj)m\ , (3-13) 

j=Q \N<\l\<2N m=l-N \m\>N I 
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where f p = 1, 4>n is an eigenvector corresponding to An, and (fj) m ( m = 0,±1,±2, ...) and (<j)N)n [n = 
0, ±1, ±2, . . .) are the Fourier coefficients of fj (j = 0, . . . ,p) and 4>n , respectively. 

4 Numerical Experiments 

In this section, wc apply Hill's method described in Section [2] to Hill's operator (|4.1|) of which the spectrum is 
exactly known, and observe that computed results are consistent with Theorems 13.41 [5751 and in Section [3J 

4.1 Hill's Operator 

As an example for numerical experiments, we consider Hill's operator [7J defined by 

s 2 = -^ + {6e 2 sn 2 (x,e)-4-e 2 }, (4.1) 

dx^ 

where sn(-,£) is the Jacobian elliptic function with modulus I (0 < I < 1). It is known [7J that the spectrum 
cr(5 , 2) is exactly given by 

<r(S a ) = k(4 -3] U K(4 0] U [a c {£), +oo), (4.2) 

where 



(a) a a {£) =£ 2 -2-2^1 -£ 2 + £ 4 , 

(b) a h (£) = -3(l-£ 2 ), (4.3) 

(c) a c (£)=£ 2 - 2 + 2y / l-£ 2 +£ 4 - 

Figure [2] illustrates a(S 2 ) for all £ with < £ < 1. 

By definition (|2.2[) of S£, p = 2 and S*^ 1 is expressed as 

with 

/i(x) = -2i/z and / (x) = {6£ 2 sn 2 (x, £) - 4 - £ 2 } + fi 2 . (4.5) 

Although the sign of the highest order derivative d 2 /dx 2 in (|4.4[) is different from that in (12.11) . we can directly 
apply the results in Section [3] to S%. Note that /o is periodic with period 2K(£), where K(£) is the complete 
elliptic integral of the first kind with modulus £. Since S% is self-adjoint and /o is analytic with d = K (vl — £ 2 ) 
in (|3.1D . Assumptions 1 and 2b are satisfied. 

For later convenience, expand the periodic term fo(x) — /i 2 in the Fourier series form 

f (x) = b + g 6, jexp + exp (~i^) } , 

where 

E(£) is the complete elliptic integral of the second kind with modulus £, and 



««>— h'-W 1 ■ 
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Fi gure 2: Com puted results of (a(S 2 ),£) for S 2 in gU (a) cr a (l) = I 2 - 2 - 
2Vl-^ 2 +^ 4 , (b) crb(^) = -3(1 - t 2 ) and (c) <t c (€) =£ 2 -2 + 2^J\-l 2 +( i . 



Then, for L = 2MK(t) with positive integer M, the eigenvalue problem S% N 4> = \(j> for 6 L2QO, L]) pcr can 
be written in the form 

N 

K{t) (Mn +J2 b ^ ((Mn-Mj + (Mn+Mj) = HMn (~N < Tl < N) (4.8) 

3=1 

with 

fyw = (v + mfffi) +bo(£) ( n = - N ---> N y ( 4 - 9 ) 

Equation (|4.8p is a matrix eigenvalue problem of 2N + 1 dimension. 

For given M and N, we can compute the eigenvalues of S£ N using Hill's method. Figure [3] shows the 
absolute errors of the computed eigenvalues corresponding to (a) in (|4.3[) for M = 2 and iV = 5, 10, . . . , 40. Here 
it should be noted that the boundaries (a), (b) and (c) in (|4.3[) are obtained as the eigenvalues of in (|4.4D 
for /i = 0. All programs used in this paper are written in C with quadruple-precision floating-point arithmetic 
and executed on a computer with the Sparc processor and the Fujitsu compiler C99. 

4.2 The Convergence Rate 

The computed eigenvalues of S® N approach to the exact ones corresponding to them with increase of N as 
shown in Figure [3] Thus the conditions (|3.3[) and p.4[) in Assumption I5al are satisfied for sufficiently large N, 
and the convergence rate (|3.9[) in Theorem 13.41 can be applied to the computed eigenvalues Ajv's. In fact, the 
exponential decay of the error with N in (|3.9[) can be found in Figure [3] Note that similar results to Figure [3] 
are obtained for the eigenvalues (b) and (c). From these, we can say that the computed results are consistent 
with (pHJl) in Theorem EHl 

For small values of £, we can show that the convergence rate (|3.9[) holds, even if the exact eigenvalues of S 2 
are unknown, as follows. For that, first, consider (|3.12j) which is the sufficient condition for (|3.3|) and (|3.4j) in 
Assumption l3al We can check (|3. 1 2[) by applying Theorem 13.31 fGcrshgorin's theorem) to the coefficient matrix 
in the left hand side of (|4.8[) . of which the diagonal elements are given by D%{£) in (|4.9p and the eigenvalues 
are real. For this matrix, Theorem 13.31 produces 

JV 

<t(S* n )c (J TO (4.10) 

n=-N 
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where 



with 



™ = { 



z G R 



^Wl<r(f)} 



r(*) 



12^ 



A'M 2 (l-g(f) 2 )(l-e^)) ! 



(4.11) 



(4.12) 



We can compute the center D%(() in (|4.9p and the radius r(£) in (|4.12p of the interval which includes 

some eigenvalues in a(S^ N )- For example, when I = 0.1 and = 0, we have 



and 



r(0.1) = 0.030- •• , 

JDg(O.l) = -3.9799 • • • , JD° x (0.1) = -2.9849 ■ 
(0.1) = 0, D% (0.1) = 4.9749- •• . 



U ±2 



Then the intervals I°(0.1) (n = 0, 1, 2, 3) are mutually disjoint and the distances between their centers are larger 
than 10 r(0.1). Moreover, from Theorem 13. 3[ we can say that I® (0.1) contains just one eigenvalue of S% N . Note 
that £ and r in (|3.12[) correspond to £>q(0.1) and r(0.1), respectively. Therefore the condition (|3.12[) is satisfied 
for the interval Iq(0.1). From these, the conditions (|3.3[) and (|3.4j) in Assumption l3"al are satisfied, and the 
convergence rate (|3.9[) holds for Iq (0.1). 

Also from Theorem 13.31 we can show that /"(0.1) and I® (0.1) contain just two eigenvalues of S£ N , re- 
spectively. Then the conditions (|3.5|) and (|3.6p in Assumption [3b] may be satisfied, and it was found that the 
numerical results are consistent with the convergence rate p. lip . For further study on this case of Assump- 
tion [3bl we have to develop the verification method for the conditions (|3.5[) and (|3.6p . 



4.3 The Error Bound 

Since the computed eigenvalues are accurate enough for sufficiently large N, as shown in Figure [3j the condi- 
tion (|3.12l) in Theorem 13.61 is satisfied for sufficiently large N. Then we can estimate the error bound (|3.13p in 
Theorem 13.61 as follows. First, £, r, L and (/o)m in p.!3p are given by, respectively, 

C = D° (£), r = r(£), L = AK(£), and (/ ) ro = (° { ™ iS ° dd) ' (4.13) 

[VLb m /2 (m is even), 

for \m\ > 1. Next, the infinite sum in (|3.13p can be bounded as 

V \(f) <2x/Z V 15 I < 24 " 2 V ^ 

|m|>JV J=W2l V ' j=\N/2~] HV ' 



(4.14) 



24^ 2 ([TV/2] (1 - gffl) + g(l)) g(l) W 2 1 
~ X(£)3/ 2 (l-<?(£) 2 )(l-g(£)) 2 

Then it follows from the inequality (|3.13p that 

^ ' \JV<|i|<2JV m=l-N 

A 247T 2 (riV/2l(l- g ffl)+g(£)) g fflW 2 ] \ 
+ + 1)^5^ (1 _ g(0 2 )(1 _ g(0)2 J - ( 4 - 15 ) 

where C, , r, L and (fo)m are given by f|4. 13|) . The computed results of the right hand side of (|4.15p are shown 
in Figure @] The absolute errors in Figure [3] arc bounded by the corresponding error bounds in Figure |U Thus 
the numerical results are consistent with Theorem [ 
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Figure 3: The absolute errors of the approximate Figure 4: The error bounds (|4.15[) of the approxi- 
eigcnvalues of cr^S® N ) corresponding to the eigen- mate eigenvalues of cr^S® N ) corresponding to the 
values (a). eigenvalues (a). 



For small values of £, we can obtain the explicit error bound (|4.15[) . even if the exact eigenvalues of S® are 
unknown. For example, when I = 0.1 and /x = 0, the interval Iq(0.1) in (|4.1ip satisfies the condition (|3.12[) 
as shown in Section T4.2I Then we can conclude that Iq(0.1) contains one exact eigenvalue of S^, and the 
corresponding error is bounded by (|4.15|) . Figure 0] shows that, when I = 0.1, the error for N > 10 is less than 
2.37 x 10" 8 . 



5 Proofs 

This section summarizes proofs of the theorems in Sections 13.21 and 13.31 Let (• , and || • H2 be the inner 
product of L2QO, L]) por and the norm induced by it, respectively. We also use || • H2 as the operator norm of 
linear operators on L2([0, L]) per with respect to the norm || • H2. 



5.1 Proof of Theorem ET3I 

To prove the theorem, we use an appropriate eigenvector </>at of the approximate operator Sr N generated by 
an eigenvector <f) of Sg. On Assumption l3"al we can employ the projection operator P\(S£ N ) to the eigenspace 
corresponding to Ayv £ a (Sp n) defined by 

Px(S^ N ) = ^i R»(S» N )dv, (5.1) 

where R v {Sp N ) = (yl — S£ j y) -1 and r\(Sg) is the directed circle with center A and radius r\ shown by 
Figured! The integral in (|5.ip is the Dunford integral. See e.g. Theorem XII. 5 in [T3]. For later use, it should 
be noted that 



Z(A(S£)) = 2^r A , (5.2) 

where l(r x (Sg)) is the length of r x (Sg), and 

\v-X\=r x , (5.3) 
rx/2< \u-X N \ <3r A /2. (5.4) 

for any v e r x (Sg). 

Then we take an eigenvector <f> with \\4>\\2 = 1 corresponding to A, and using <f> and the projection P\ N {S^ N ) 
for X N G a(S£ N ), we can obtain an eigenvector 4>n with 1 1 jv 1 1 2 — 1 corresponding to Ajy- The feasibility of 
such procedure is guaranteed by the following lemma, which is based on the idea in [T] . 
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Lemma 5.1. Let satisfy Assumvtions\l\ and \2al A G o~{S£) and Ajy £ °~{Sp n) satisfy Assumption ] Sal Then 
for an eigenvector <fi corresponding to A we have 



U-PxAs;, N )<t>h < -W$ - s^th- 

Furthermore, if \\4>\\2 = 1 and there exists a positive integer No such that for any N > iVo 

2 



'A 



\\(S^-S^\\ 2 <1, 



(5.5) 



(5.6) 



we have P Xn (S^ N )<b ^ for N > N . 

Proof. Under Assumption 2a, we can use the contour 7a (S^) to define P\(S£) and P\ N (Sp N ). By the residue 
theorem, we have P\(S£)4> = <fi. Then it follows from (|5.2[) that 

H-PxAS^h = 11^(^)0-^(^)0112 

i 



2?ri 



{R v {s>;)-R v {s; tN ))^Av 



< r x max \\(R V (S£) - R v (S%, N ))<f>\\ 



(5.7) 



As for \\(R u (Sg) - Rv(S»jf))<l>\\2 in (& we have 



||(iZ„0S£) - R v (S£ N ))<f>h = II - - SfaHvI-Sp-^h 

= II - Sfay^S* - S£ n )(p A)-V|| 2 

< |i/ - Ar 1 ||(„7 - S^rla ||(S£ - 5 p %)0|| 2 
1-1 



P,J 



(5.8) 



= \u - X\- 1 dfctfa a^))- 1 || (S£ - S£„M| a , 
where dist(i>, o-(S^ N )) is the distance between ^ e C and c(S^ N ) defined as 

dist(i/,<7(5^))= inf >-£|. 

f e<T ( S p,Jv) 

The last equality in (|5.8|l follows from the self-adjointness of S£ N . Furthermore, by (|5 .3|) and (|5.4[) we have 

|v-A|=r A , (5.9) 

for any v £ 7\(S^) and 

min dist(v, a(S% N )) = mm \v-\ N \>r\/2. (5.10) 
Combining (jS~5)l . ([53)1 and ([5710)1 we have 



max \\(RAs^) - RAs» N ))<i>h < 11(3? ~ SZnWU- 



(5.11) 



Then we can deduce the conclusion (|5.5[) from (|5.7[) and (|5 . 11 [) . 

Finally, if ||0|| 2 = 1 and ([575]) holds for N > N , it follows from ([575]! that ||0 - P Ajv (S^^jb < 1. Hence 



P\s (S^jv)^^O holds true. 



□ 



Here we present a key lemma for the estimate of the convergence rate. Part of this lemma is also based on 
the idea in [T|. 
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Lemma 5.2. Let S£ satisfy Assumptions]]] and \2al X 6 cr(S^) and Xn S a (Sp n) sa -ti s fy Assumption \Sal 
Furthermore, for an eigenvector corresponding to X with \\<p\\2 = 1? assume that there exists a positive integer 
No such that for any N > Nq 



-H(S£-S£ w M| 2 <i. 

Then for N > Nq we have 

\x-x N \< (3+^) m - s» N )ct>\\ 2 . 

Proof. By Lemma T5. 11 for N > Nq we can define 



\P\ N {Sp 



Then is an eigenvector corresponding to Ajv with || || 2 = 1. Then noting A = (S£(, 
( S p,n4>n, 4*n)2 we have 

A - X N = (S^0, 0) 2 ~ (Sp <N (f>N, 0^)2 

= (S£(j>, - 0jv) 2 + ((S£ - S£ N )<j), <t>Nh + (<t> - <Pn, S^ N (f> N ) 2 , 

where the last equality follows from the self-adjointness of S£ N . Therefore we have 



IA-AjvI < (\\\ + \\ N \)\\<i>-<l>Nh + KS£-S» N ) 



Here ||0 



'■>N\2 



is bounded by 



110- fah < H-Px N (s^ N )4>h + WPxAS^-M 



Since \\Px N (S£ 



N) 



1 2 in (|5.17p is bounded as 



1 \\PxAS%^h) PxAS ^ 



we have 



< 11^(^)0-0112, 



||0-0iv|| 2 <2||0-P Ajv (5^ JV )0|| 2 . 



Applying the inequality (|5.5|l of Lemma 15.11 to (|5.18[) , we have 

110- <Mh< ^11(^-^)0112. 

Then combining (|5.16[) and f)5 . 19[) . we have 

IA-AjvI < rA(|A| + |Aiv|) + l) ||(^-^)0|| 2 . 
Finally, noting |Ajv| < |A — A jv | + | A | < r\/2 + |A|, we have the conclusion. 



(5.12) 
(5.13) 

(5.14) 

) 2 and A jv = 

(5.15) 

(5.16) 
(5.17) 



(5.18) 



(5.19) 



□ 



Remark 5.1. Assumption]]^ i.e. the self-adjointness of is used (i) to estimate the norm of the resolvent 
of the approximate operator as \\(yl — •S'p jv ) _ 1 1 1 2 = dist(i/, <r{S^ N )) _1 in (|5 . 8[) . (ii) to make an approximate 
eigenvector as (|5.14[) and (Hi) to deduce the equality (|5.15l) . 
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What remains is to estimate \\(S^ — S£ N )<fr\\2 in the RHS of f)5.13[) . Since 

II - s» N )<t>h = US'; - PnS^ + (P n s; - P N s;P N )q>\\ 2 

< ||(I - P N )S^\\ 2 + \\PnS$(I - Pn)<P\\2 

< (|A| + \\P N S^h)\\(I - P N )4,\\ 2 , 

we consider estimation of \\PnS£ 1 1 2 and || (I— Pn)4>\\2- Their estimates are given by the following lemmas, whose 
proofs are shown in Appendix El It suffices to consider the case of Assumption l2"al for ||P/vS^||2, whereas the 
both cases of Assumptions l2"al and I2bl need to be considered for |j(I — Pn)4>\\2 which depends on the smoothness 

of (f>. 

Lemma 5.3. Let S£ satisfy Assumvtion \2al Then there exits a positive real constant C depending only on p 
and fo, fi, ■ ■ ■ , fp-i such that 

\\PnS;\\2<CNP+^ 2 . 

Lemma 5.4. Let S£ satisfy Assumption \2a\ and (f> £ H p ([0, L]) pcl - be an eigenvector of S p l corresponding to 
A G c(S^) with ||</>||2 = 1- Then for any positive integer q there exists a positive constant C depending only on 
4> and q such that 

Kl-PN)<l>h<C(l + N)-«. 

Lemma 5.5. Let S£ satisfy Assumption{2j^for d > and <fi € H p ([0, L]) pel be an eigenvector of correspond- 
ing to A £ cr(5'p) with \\4>\\2 = 1- Then for any e with < e < d there exists a positive constant C depending 
only on <f>, d and e such that 

II (I - pN^h < C cxp (- Md L ~ £) N\ , 

where L is the period of the coefficient functions of S£ . 

Here we prove Theorem 13.41 
Proof of Theorem \S.4\ Since it follows from Lemmas 15.31 and 15.41 that 

\\(S^ — S; N )4>\\2 — » as N — > 00, the condition ()5.12[) in Lemma I5~2l is satisfied. Therefore we can deduce the 
conclusion (|3.8[) from Lemmas 15.21 15.31 and 15.41 and p.9[) from Lemmas 15.21 15.31 and 15.51 □ 



5.2 Proof of Theorem 1331 

We use almost the same methods as the ones of Lemmas l5.ll and l5.21 On Assumption l3bl we consider Pa n {S^ N )(f> 
for an eigenvector tj) corresponding to A € v{Sp N ), where An = {Ajv,i, . . . , Xnm}, 

PaAS^.n) = / R^S^dft, (5.20) 

and T\(S^) is a directed circle with center A and radius r\, which contains A^. Here it should be noted that 
(|5.2p . (|5.3[) and (|5.4[) can be used for Assumption [3b] when \jy is replaced by Aiv.i {i = 1,2, ... ,k). In this 
case PA N (Sp N )4> is not an eigenvector corresponding to a single eigenvalue of S£ N but a linear combination of 
eigenvectors corresponding to Aat.i, . . . , Ajv^. Then we can estimate the convergence rate of one of them after 
appropriate rearrangement of them explained later. 

Lemma 5.6. Let satisfy As sumptions [7] and \2a[ A G c(S^ ) and \n,i, • ■ • , ^N.k € a (Sp n) sa ti s fy Assump- 
tion WR Then for an eigenvector tfi corresponding to A we have 

H - Pa n (S^ n )<P\\2 < ^-\\(S£ - S» N )<l>h, (5.21) 
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where An = {Ajv,i, • • • , ^N.k}- Furthermore, if \\4>\\2 = 1 oxid there exists a positive integer Nq such that for any 
N > N 

^\m-S^ N )4>\\ 2 <l, (5.22) 

we have P An (S^ N )(f> ^OforN>N . 

Proof. Noting the contour r\(S£) = dB\{r\) surrounds A and Ajv,i, • • • , Ajv.fe in the case of Assumption I3b[ we 
can use r\(S£) to define P\(S£) and Pa n (S^ n ). Then we can obtain the conclusions in a similar manner to 
Lemma 15.11 □ 

Lemma 5.7. Let S£ satisfy Assumptions^ and \2al A <E cr(S^) and Ajv,i, • ■ • , ^nm G a {Sp n) sa ti s fy Assump- 
tion\3b[ Furthermore, for an eigenvector (j> corresponding to A with \\4>\\2 = 1, assume that there exists a positive 
integer Nq such that for any N > Nq 

^ll(S£-S£ w M|a<l. (5.23) 
Then for N > Nq there exist a choice of an index i^ € {1, . . . , fc} and an integer k' N with 2 < k' N < k such that 

|A - \s,is\ < [Kx,N \\{S% - S^ N )<Ph) , (5.24) 

where 

k 2 max jl, (|A| + ^) 2 "| max {l, \\S^ N \\l) 



K X , N = (4|A|+r A ) 

max] 1,(11^112 + ^ 



2 



max < 1 , 



2 



In particular, if there exists N\ > A^o such that K\.n \\ {Sp — N )4>\\2 < 1 for N > Ni, the integer k' N in (|5.24D 
can be replaced by k. 

Proof. If there exist i,j £ {1, ...,k} with i ^ j such that \N,i = Ajvj, we reserve only one of them, and 
consequently we obtain the set of indices O^r such that {Xna | i <£ ^n} is identical to {Aa^i, . . . , \n,Ic} as a set 
and the elements of it are mutually distinct. Then we renumber the elements of {Xn,i \ i € and let it be 
denoted by {\n,i, ■ ■ ■ i Ajv.fc' }, where k' N = #^at. Without loss of generality we can assume that 2 < k' N . 

It follows from Lemma 15.61 that Pa n {S^ n )4> is a nonzero linear combination of eigenvectors </>jv,i, ■ • ■ , 't'N.k^ 
corresponding to Xn.i, ■ ■ ■ , A n ,k' N for N > Nq. Then there exists a sequence of coefficients cln,i, ■ • ■ , o,N,k' N <= C 
with (ajv,i, • ■ • ,a>N,k' N ) (0, ■ • ■ , 0) such that 

PA N {Sp N )4> = ajv,ife,i H h a N ^ N 4> N M N (5.25) 

for N > Nq. Here we define ijv(l) as a minimizer of |A — Ajv,i| for i with ajv,j 7^ in (|5.25|) . i.e. 

|A - \ N ,i N (i)\ = min{|A - Ajv.il | i : ajv,* 7^ in dE23}. (5.26) 

Moreover, we define {^(2), . . . , i^(k' N )} as a permutation of {1, . . . , k' N } \ {^(1)}. Consequently, we obtain 
the ordered set {Ajv,jjv(i), • • • , ^N,i N (k' )} for N > Nq. For simplicity, let Xnj, 4>nj and oj\r,j denote X}f,i N (j), 
4>N,i N (j) and a>N,i N (j) f° r j : = 1, • ■ • , k' N , respectively. In the following we consider (|5.25|) under this renumbering. 

Here we begin the estimate of the convergence rate of Xn,i- Applying the operator (St 1 N — Xn^I) • • ■ (S£ N — 
Ajv.fc' -0 t° the both side of (|5.25p . we have 

9n(S£ >n ) Pa n {Sp tN ) <j> = a Ni i ffjv(A^,i) <j>N,i, 
where gn is the polynomial defined as 

9n{x) = (x- Ajv, 2 ) ■ • • (x - Ajv.fc' )• (5.27) 
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Because ajv,i ^ and \n,i 7^ Ajv,?: (i = 2, . . . , fc^v), we can define a normalized eigenvector 0jv corresponding to 
Ajv,i as 

<t>N = C^ 1 <?;v(S p %) P An (S» n ) </>, (5.28) 

where 

C N = sign(. 9Ar (A)) || 3jv(^ iJV ) PaAS'n) (5.29) 
Using in (|5.28[) . in the same manner as (|5.16p . we can derive an estimate: 

|A - X NA \ < (|A| + \Xn,i\)H- <t> N h + \\(S$ - S» N )4>\\2. (5.30) 
Multiplying the both side of (|5.30p by |<?jv(A)|, we have 

I] l A " V*l < (|A| + |Ajv,i|) \9n(X)\ U - cf> N \\ 2 + |. 9A r(A)| ||(S£ - S^\\ 2 . (5.31) 

Noting Assumption [3b] and (|5.26p , we can deduce from (|5.31[) that 

|A - A Nil | fe « < (|A| + IAjv.iI) \g N (X)\ H- fah + (rxm^US^ - S» N )<t>\\ 2 . (5.32) 
Then what remains is to estimate |<?at(A)| \\<f> — <$>n\i- By the triangle inequality, we have 
\g N {X)\ U-MV< l|5Jv(A)0 - C N (f> N \\ 2 + \\(Cn - 9n(X)) <t> N \\ 2 , 

and 

||fl r Jv(A)0 - C N 4> N \\ 2 

= \\g N (S£) Px(S£) <f> - g N (S^ N ) P An (S» n ) 0|| 2 

< \\{9n(S^ - gN(S* N )}<f>\\ 2 + \\9n(S^ n ){Px(S^ - Pa n (S^ n )}^\\2. (5.33) 

Therefore setting 

Ei,n = || (Cjv - 9n(X)) 0jv H2, 

^2,iV = ||{3Jv(^)- 5J v(^ JV )}0||2 1 

we haw 

| 5A r(A)| 110- <f> N \\ a < E hN + E 2 . N + E 3 . N . (5.34) 

We can estimate E\,n, -E 2 ,jv and S3,jv as follows. Using (|5.28[) and (|5.29p . we have 

E hN = \C N - g N {X)\ = \C N - sign(g N (X))\\g N (X)(j)\\ 2 \ 
= \\\C N (j) N \\ 2 - \\gN(X)4>\\ 2 \ < \\C N 4> N - g N (X)4>\\ 2 

< E 2 , N + E 3 , N , (5.35) 



where the last inequality is due to (|5.33[) . Noting that <j) is an eigenvector of S£ corresponding to A, we have 
MSg) - g N (S^ N )}^ = {g N (X) - g N (S» „)}<!> 



few 1 



= ]T a(z,A Ar )(A fc «- 1 -V-(^ JV )' £ «- 1 - i ) 

i=0 

= <x(i,A N )/3(i,X,S» N )(XI-S» N )<f> 

k N —2 

= J2 A ~) P& A, S^ N ) (S£ - S» N ) 0, (5.36) 
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where a(i,A N ) = (-l) , E 2 < 3l <-<3,<fc 
fore noting 



X NJl ■ ■ ■ X NJi and (3(i, A, S») = E^ 2 * A*^- 2 "*- 1 ^ „)». There- 



a(i,A N )\< |A| + -£ <max 1, |A| + 



AT 2 



(i, A, 



,JV 112 



^)|| 2 <max{l,|A| fc '«} Y, \Ki 

1=0 

{l,|A| fc «} k' N max {l,||^J^} 



< max • 



and f|5.36[) . we have 

£ 2 ,w < (fc^) 2 max<| 1 
Using (|5.21[) in Lemma |5~5I we have 



max 



(5.37) 



^ < A \\s* N - \ N ,ii\\ 2 \\{Px(sp - Pa n (S^ n )M2 



ii=2 



< 



jQK N h + |Aiv,i|) 110-^(5^)0112 



i-t=2 



< l|5 p ^|| 2 + |A| + f -11(^-5; 



Finally, noting k' N < k and combining (|5l?2"j) . (|5TM|) . ([535]) , (|5T57|) and (j535|) . we have 
|A- A^il^ 

< 2(|A| + IA^xD^jv + E 3iN ) + (rxfif"- 1 ^ - S» N )<f>\\ 3 . 



< 



k'^m&xl 1, ( |A| + — 



(4|A|+r A ) 



-Jl m ax<jl,(||5; iJV || 2 + |A| + ^ 



2 k 



( y 



I6S£-S£ 



p,N) 



Thus we obtain the conclusion. 
Here we prove Theorem 13.51 



(5.38) 



□ 



Proof of Theorem \3.5\ In a similar manner to the proof of Theorem 13.41 we can deduce the conclusion (|3.10D 
from Lemmas lCT Pl andlBTTl and (|3~TT|) from Lemmas lQISn andlSTfl Note that ||S£ 



< WPNSHhWP, 



2 -nv 2 



l^|| 2 . 



□ 



5.3 Proof of Theorem [3761 

To prove the theorem, we use a similar technique to Lemmas 15.11 and 15.21 in which we exchange the roles of 
4> and 4>n- That is, we take an eigenvector (f>N with ||0j\r||2 = 1 corresponding to An, and using (f>N and the 
projection P\(Sg) for A G c{Sp), we obtain an eigenvector <f> with ||0j|2 = 1 corresponding to A. First, we 
present a lemma guaranteeing the feasibility of such procedure, which corresponds to Lemma 15.11 We omit its 
proof since it is proved in almost the same manner as Lemma 15.11 



Lemma 5.8. Let satisfy Assunwtions\l\ and \2a[ A € o-(Sg) and Ajv € cr{Sp N ) satisfy Assumvtion\3a\ Then 
for an eigenvector <pN corresponding to Ajv we have 



\<Pn - Px(S^ N \\ 2 < -||(££ - S^Nh 



(5.39) 
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Furthermore, if \\4>n\\2 = 1 ond 



^II(S£-S P %)<M2<1, (5.40) 



we have P X {S$)4> N ^ 0. 

Next, we present a counterpart of Lemma 15.21 whose proof is also omitted. 



Lemma 5.9. Let S£ satisfy Assumptions]]} and \2a[ A € o~(S£) and Xn £ a(S p N ) satisfy Assumvtion \3al 
Furthermore, for an eigenvector 4>n corresponding to Xn with \\<Pn\\2 — 1; assume that 

^ll(^-^)^||a<l. (5.41) 



Then we have 



Here we prove Theorem 13.61 



|A - A,v| < (^3 + || ( 5 m _ S» N )d> N \\ 2 . (5.42) 



Proof of Theorem \3.b\ It follows from the assumption (|3.12D . Theorems 13.11 and 13.21 that the set of the accumu- 
lation points of the sequence {Xn} coincides with cr(S£) H B^(r). Moreover, the elements of <j(S£) \ B^(r) do 
not exist in the interior of B^(9r). In the following we set A^ >r = cr(S£) H B^(r). 

First, we show that {Xn} is convergent and A G A^ iT is uniquely determined as the limit of {A^}. To prove 
{Xn} is Cauchy, we choose an arbitrary integer M > and estimate |Ajv+m — Ajv| for sufficiently large N. For 
an eigenvalue A £ A^, r , we take an eigenvector <j> with ||0||2 = 1 corresponding to A. Then <f> satisfies 

P A( JS^=^-£ (vI-Sg)- 1 <t>dv=^- i £ {v-X)- 1 4>dv = <t>, (5.43) 

where J^. r is the boundary of B^(2r) with counterclockwise direction. Note that (|5.43|) holds true regardless of 
the number of the elements in A^ <r . Using this 4>, we define eigenvectors corresponding to Xn+m and Ajv as 

<pN+M = Tj-5 77^ r7]r- P A« + Af(' S 'p,JV+A/)^ 

\\"Xli+u\. :} p,N+M) ( P\\2 

and 



>N = 



-PxAs, 



p,N) 



respectively. Here we use again r^ r to define P\ N+M (S p N+M ) and Px N (S p N ). The availability of the eigen- 
vectors can be guaranteed by a similar manner to Lemma |5. II Then noting S p N = PnS^Pn etc. we have 

Xn+m - Aw = (S ptN+M </>N+M,<j>N+M) — {S PjN+m ^n,4'n) 

= {S^n+M^N+mAn+M - 4>n) + {S p!N+M ((f)N+M - 4>n),4 i n) 
= (S£ N+M <f>N+M,<f>N+M ~ 4>n) + (<f>N+M - 4>N, S p N+M 4>N) 

and therefore 

|Aat +A / - Ajv] < 2\\S^ N+M UU N+M - <f>\\ 2 + \\<f> - fah). (5.44) 
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Using a similar manner to Lemma 15.21 Lemma 15.31 and Lemma 15.41 we can show that the RHS of (|5.44p is 
0(N~ q ) for any positive integer q. In fact, \\S£ N+M \\ 2 = 0(iV p+1 / 2 ) by Lemma IQ1 and 

110 - 0isr|| 2 < 2 PxAS^h = 2 ll^c.OW ~ PxJSp.N^h 
< 4r max \\{yl - S^ N )-\S^ N - S$)(yl - Stf-^h 

= 4r max \\(yl - S» N )-\S» N - 5£)(i/ - A)"V||2 



< 4r max 



ver C:T \v — A| 1 — A 7v | 
<^II(^-W| a . 



Thus we can show that |Ajv+m — Xn\ — ► as N — > oo and A € ylf >r is uniquely determined. Combining this 
fact and |A — Ajv| < 2r, we can show that Assumption I3al is satisfied for {An}, A, and r\ = 4r. In fact, for any 
A G a(Sg) U £7(5^) \ {A, Ajv} we have 

|A-A|> |A-C|-|A-C| >2r A . 
Next, we show the error bound (|3.13p . By Lemma \5. 91 it follows from (|5.42j) and |A| < r + |£| that 

3 ICI\. nrau cm w n (5 45) 



\\-X N \< 



for A £ <?(S£) H B^(r). Furthermore, we have 

(S£ - )^ = (/ - P N )S^ N 



p l+N 



= E \ E E {fi)m{4>N)l-m^7f ( -i 
\1\>N \j=0m=l-N V L 



1 / . 2tt(/ - m) 



.2irl 

cxp | — 1— —X 



and therefore 



p i+w ^ 1 / 2ir(l — m) x 

E E (fj)m($N)l-m-^= (-i 



j=0 m=l-N 
p P 



L 



^4f(¥T££ E |(/i)m(to-™ 



VL \ L J 



j=0 \1\>N m=l-N 



Considering the case \l\ > 2N, we have 



l+N l+N 

E E |(/i)m(0iv)i- m |< E E |(/>)m| <(2JV+1) E |(/i)r 
|Z|>2iVm=i-iV |;|>2ATm=;-Af |m|>JV 

Combining (|5.45l) . (|5.46p and (|5.47p . we obtain the error bound (|3.13|) . 



(5.46) 



(5.47) 
□ 



6 Concluding Remarks 

We have considered the convergence rates and the explicit error bounds of Hill's method, which is a numerical 
method for computing the spectra of self-adjoint differential operators S p in with periodic coefficient 

functions fj. On the assumption (Assumptions [3"al and l3b|) that the computed eigenvalue A at is close to the 
exact one A, it is shown in Theorems 13.41 and 13.51 that the convergence rate of the computed eigenvalue is all 
order polynomial in the case of fj being C°° as shown in (|3.8[) and (|3.10[) . and exponential in the case of fj 
being analytic as shown in (|3.9|) and (|3.11[) . In addition, even if the exact eigenvalue A is unknown, it is shown 
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in Theorems 13.41 F5T51 and I3~6l that, if the condition (|3.12[) is satisfied, we can obtain the convergence rate using 
p.8p - (|3.11[) and the explicit error bound using (|3.13[) . There are some cases in which the condition (|3.12[) can 
be checked using Gershgorin's theorem. These theorems are proved using the Dunford integrals P\(S£ N ) in 
()5.1[) and P\(S£) in the proof of Theorem 13.61 which project an eigenvector to the corresponding eigenspace. 
This integral is suitable for proving these theorems, because the boundary of the neighborhood of an eigenvalue 
can be directly used as the contour of the integral. Numerical examples using Hill's operator (|4.1[) support these 
theoretical results. 

As described in Remark 15.11 self-adjointness of the operator enables us to estimate the norm of the 
resolvent as shown in (|5.8[) . and to project an eigenvector to the corresponding eigenspace using the Dunford 
integral. For the case of non-self-adjoint operators, we have developed some other approaches, which will be 
reported somewhere else soon. 
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A Proofs of Lemmas 15. 3L 15. 4L and 15.51 



First, we prove Lemma [ 

Proof of Lemma COI We set f p = l for simplicity. Let e m (m = 0, ±1, ±2, . . .) be the Fourier bases defined as 



6m (^O 



cxp — 1 



. 2irm 



L 



and let (fj) m (j = 0,1,..., p, m = 0, ±1,±2, ...) denote the Fourier coefficients of fj (j = 0,1,..., p), 
i.e. (fj)m = {fj, e m )2- Then it is a standard fact that under Assumptionl2"aI for j = 0, 1, . . . ,p, m = 0, ±1, ±2, . . . 
and any positive integer q, there exists a positive real constant Cj tQ depending only on fj and q such that 



\(fj)m\ < Cj, q (1 + \m\)- q . 



We estimate the norm of 



N p ( oo 
n= — Nj=0 \m= — oo 



. 27r(rt — m) 



1pn—m ( 6n(*^) 



for JV > 1 and a nonzero vector ^ € £2([0, £]) pcr . Setting 



Cj(n,m) = (fj) m - 



. 2it(n — m) V 



and q = p + 1 in (|A.1[) . for j = 0, 1, . . . ,p, m = 0, ±1, ±2, . . . and n with |n| < JV, we have 



) 



2tt 



|n| + |m| 
JV 



< Mr <W (i + H)^ 1 i + ^ 



JV 



and therefore 



Combining (|A.2[) and (|A.3|) we have 



|cj(n,m)| < 



27r\ Cj, p+ i 

Ty i + |m| 



JV J . 



ip^ii 2 < E 



JV 



E E Cj{n,m)^ n - m ) e 

-TV \m— — oo 
oo 

Cj(n,m)4> n - 



P N 

= E E 

j=0 I n= — TV m= — oo 
p ( N / oo 

<E E E i c i(«.-)i 2 

j— v n— — N \m— — oo / \m — 

s£ £ (|) <S*.W E 

3=0 ln=— iV v / \m=-( 

<EQ, P+1 (2JV+1) 1 / 2 JV^||^|| 2 , 

3=0 



1/2 



E l^n- 



(1 + H) 2 



1/2 



HV'lli 



1/2 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
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where Cj.p+i (j — 0, 1, . . . ,p) are bounded constants defined as 

1/2 



c 3 , P+ i - c jtP+1 I e | j 

v 7 \m— — oo v / 



Hence it follows from (|A.4|) that 

||iV?£|| 2 < E Q, P+1 (2A + l)V2j\P. (A.5) 
Thus we have the conclusion. □ 



Next, we prove Lemma 15.41 and 15.51 with the aid of the following two lemmas. Since these lemmas are 
fundamental in the theory of differential equations, we only give the outlines of their proofs in Appendix [Bl 

Lemma A.l. Let S£ satisfy Assumvtion \2a\ and <j) G H p ([0, L]) por be an eigenvector of S£ corresponding to 
A G o-(Sg). Then cj> is C°° . 

Lemma A. 2. Let S£ satisfy Assumption\2j^ for d > and cf> G H p ([0, L]) pcl be an eigenvector of S£ corre- 
sponding to A G a(S£). Then <f> can be uniquely extended to an analytic function on T>d with <p(z + L) = 4>{z) 
for any zeC, where T>d is defined as (|3.1j) and L is the period of the coefficient functions of S£. 

Proof of Lemma \5.4\ Since <f> is C°° by Lemma lA.U in a similar maimer to (|A.1[) . we can use the fundamental 
fact for the Fourier coefficients that there exists a positive constant C^^ depending only on (f> and q such that 

\k\ < cv, 9 (1 + Inl)" 9 " 1 

for any integer n. Then we have 

11(7-^)^111= E \^\ 2 < E Cl q (l + \n\r 2 "- 2 

\n\>N |n|>JV 

C2„ _ i cl 



< <t>,Q 



^ (i + M) 2 - (i + n) 2 i ^ ri + |j 



(1 + N y q ^ (1 + | n | )2 - (1 + N?q ^ (1 + | n | )2 ■ 
Thus we obtain the conclusion. □ 



Proof of Lemma \5.5\ Since 4> can be regarded as an analytic function on T>d with <f>(z+L) = <f>(z) by Lemma lA.2[ 
there exists a positive constant C^^.e depending only on <j) 1 d and e such that 

\k\ <C Me exp(-^-\n\j (A.6) 

for any integer n, where d e = d — e. In fact, we can show this estimate by the Paley- Wiener type argument. 
Noting the periodicity <fi(z + L) = <f)(z), by Cauchy's integral theorem we have 

1 f L , ( .27m 



VI 



J 4>{x) cxp ^— i — x^j dx 



1 f L , , ( .2nn 



(a; — sign(n) d e i) exp I — i— — {x — sign(n) d e \) I dx 
vL Jo \ L ) 

( 2nd, A 1 f L . n ... / .27rn \ 

= exp I — -|n|l — = y (a; — sign(n) a e i) exp I — i— — x I dx, 

where sign(n) is the sign of n. This expression implies (|A.6[) . Therefore we have 

. „ _ „ 2C? rf ,exp(-2d £L (A + l)) 

|n|>JV |n|>JV Pl £ ' Lj 

where d 6i i = 2nd £ /L. Thus wc obtain the conclusion. □ 
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B Proofs of Lemmas IA.1I and IA.2 



Proof of Lemma \A.l[ Since <f> is an eigenvector of S£ corresponding to A, we have 



p-i 



3=0 



Then it follows from (f> G H p ([0, L]) pel that ^ G i^QO, L]) per . Hence G H p+1 {[0, L]) pcr . Iterating a similar 
argument, we have <p G ^ 9 ([0, ^])per for any positive integer q. Then the conclusion follows from the fundamental 
embedding relation: H q ([0, L]) pei °->- C 9_1 ([0, L]) per for any positive integer g, where C" 3_1 ([0, L]) pGr is the set of 
periodic functions on [0, L] which arc g — 1 times continuously differentiablc (see e.g. [3J Proposition 7.5.4]). □ 



Proof of Lemma \A.2[ We choose an arbitrary e with < e < d and set d e = d — e. Then /o, /i, . • ■ , fp—i are 
analytic on T>d e and continuous on T>d E - Hence there exists a positive constant M £ such that |/j(;z)| < M e for 
all j = 0, 1, . . . ,p — 1 and z G T>d e ■ 

By Lemma fA. 11 we can assume the existence of an eigenvector <fi G C°°([0, £]) pe r corresponding to A G <r(S£). 
Then for the fixed A and each xq G R, we consider the equation S^ip — XtJj as an initial value problem with the 
condition ip^ n \xo) = (j)^ (xq) (n = 0, . . . ,p — 1). By the matrix F(z) defined as 



F(z) = 



/-f P -i(z) -fp- 2 (z) 







V o 

and the vector w = (0( p_1 ' (x ), <^ p_2 ) (x ), 



-A(z) -(f (z)-X)\ 




1 / 

(xq)) t , the initial value problem above is expressed as 
w'(z) = F(z)w(z), w(xq) = wq, (B-1) 
where w(z) corresponds to (z), i/j( p_2 )(z), . . . , -0(z)) T . The equation (jB.ll) is equivalcntly expressed as 

ti>(z) = w + [ F(s)w(s)ds. (B.2) 



By the fundamental method of successive approximations Theorem 2.3.10, it is shown that the equa- 
tion (|B.1[) has a unique analytic solution w on T>d e for each xo G R. Then what remains is to prove the 
consistency of the solutions for different xo's and the periodicity of the solution. First, note that for any Xq the 
solution w is identical to v — (d^ p ~ 1 \ d^ p ~ 2 ^ 6) T on R. This follows from the condition 



\w{x) — v(x) 



F(s)(w{s) -v(s))di 



< sup ||F(z)|| r 

z€T> dE Jx 



where || • |j is a norm of R p . Then if w\ and W2 are the solutions for x^ and x 1 ^ 1 , respectively, we have 



w(s) - v(s)\\ \ds\, 



,.< 2 ) 



(B.3) 



Wi(z) — W2{z) ~ wi(x^) — w 2 {x\^ ! ) + I F(s)wi(s)ds 



(2) 



i ( ( 2 h 



(1) 



X W 
(2) 



(2) 



F(s)w2(s)ds 



F(s)w 2 (s)ds+ / F(s)(wi{s) - w 2 {s))ds 



= v(x { Q ) )-v(x i o ) )- I a F(s)v(s)ds+ I F(s){ Wl {s) - w 2 {s))ds 

c x 



F(s)(wi(s) - w 2 (s))ds 



1 In [9| Theorem 2.3.1], the general form w' = F(z,w) is considered. Here we consider the special case that F is linear with 
respect to w. 
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for any z £ T>d e , which implies w\ = W2 on T>d € by the same criterion as (|B.3[I . Next, the expressions 

/•Z + L rZ 

w(z + L) = v(x + L) + F(s)w(s)ds = v(x ) + F(t + L)w(t + L)dt 

Jxq-\-L Jxq 



= v{x )+ / F(t)w(t + L)dt, 

Jx 

w(z) = v(x ) + [ F(t)w(t)dt. 



guarantees the periodicity of w also by the same criterion as (|B.3I) . 

Finally, since the all arguments above hold true for any e with < e < d, we obtain the conclusion. 
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